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' Abstract 

The final state interactions in tt^tt^tt'' decays are considered 

using the methods of non-relativistic quantum mechanics. We show how 
^ \ to talce into account the largest electromagnetic effect in the analysis of ex- 

■ perimental data, using the amplitudes calculated earlier. We propose the 

relevant expressions for amplitude corrections valid both above and below 
the two charged pion production threshold M^o^o = 2m^±, including the 
average effect for the threshold bin. These formulae can be used in the 
procedure of pion scattering lengths measurement from M^o^a spectrum. 

O ; The experiments DIRAC and NA48/2 at CERN SPS are able to mea- 

r~| ' sure the tttt scattering lengths difference ao — 02. Its value is predicted by Chiral 

CLii Perturbation Theory (ChPT) with a high accuracy [3] (ao - 02 = 0.265 ± 0.004 

in units of inverse pion mass.) Thus the extraction of ao,a2 from experimen- 
I tal data with comparable precision become an important task. Recently the 

• NA48/2 experiment has discovered an anomaly (cusp) in the invariant mass 

M^o^o spectrum of 7r° pairs from — > Tr^Tr^Tr" decays. This anomaly is posi- 
, - tioned at the charged pions production threshold Mc = 2m^± , and has been ex- 

rN I plained by N.Cabibbo [4] as a result of charge exchange reaction tt+tt" n'^n'^ 

^ . in final state. It turns out that investigation of these decays below the charged 

pions threshold Mc allows one to extract the scattering lengths with a precision 
not accessible in other current experiments. 

Further development of the model has been done in the works [S1E1I7], where 
the pions final state interactions were considered with the second order precision 
in scattering lengths terms. Approach developed in [S] provides the possibility 
to extract ao,a2 from experimental data [2] with a precision compatible to 
theoretical prediction one. Nevertheless there are two problems, requiring a 
special consideration, that can be crucial for scattering lengths extraction from 
experimental data in the vicinity of the threshold. 

The first problem is the accounting for a higher orders in scattering lengths 
terms and their impact on the value of ap, 02. The estimation made in [S] leads 
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to the corresponding 5% error for ag — 02. In comparison with the precision 
of ChPT prediction and with the experimental statistical error this is a rather 
noticeable uncertainty, and its decreasing is very desirable 0- 

Another effect that has an impact on ao,a2 extraction from experimental 
data is the electromagnetic interaction of pions in the final state. 

This issue is not a trivial task [TI] [T^], as Coulomb interaction between 
charged pion^ below threshold Mc leads to formation of bound tt+tt" states 
(pionic atoms A27r±), and construction of them in the framework of perturbation 
expansion is a doubtful issue. 

To solve this problem we use the technique of non-relativistic quantum me- 
chanicH. Leaving the strict derivation for elsewhere, we cite here the main 
result for amplitude of the decay under consideration 

T (1 + it;o/o)ro + 2iw/j;T+, fa = aoo/D; f.^ ^ a^/D, 

D — {1 — ivoaoo){l ~ 2iva-i ) + 2vova^ (1) 



Here mo,m are the Tr",n^ masses, and vq = - — — -;v — - — ^ are 

the neutral and charged pions velocities respectively. and To are the "un- 
perturbed" matrix elements of tt^tt+tt^ and — ^ Tr^Tr^Tr" decays. 

The inelastic and elastic aoo, an pion-pion scattering amplitudes in the iso- 

topic symmetry limit are E] a^, = (oq — a2)/3;aoo = (oq + 2a2)/3;a_| = 

(2ao + 02)76. The replacement — *■ fi has small numerical impact on results 
of previous calculations done according [SJ [5], but is crucial for inclusion of 
the electromagnetic interactions under threshold, where the formation of bound 
states (tt+tt" atoms) take place. 

The expression ([TJ includes all successive elastic and inelastic interactions 
in the tttt system (to all orders in scattering lengths). Their interaction with a 
spectator pion can be taken into account at two- loop level O [5] . 

As it was discussed earlier jTU], to include the Coulomb interactions in the 
frameworks of considered approach, it is enough to make the simple replacement 
in H]): 

cx 

V T = iv — a [logf— 2TOmro) -f 27 -f — ; ^ = — (2) 

2v 

where 7 = 0.5772, a = 1/137 are Euler and fine structure constants, whereas 
— '^'"rfg"^^'* is digamma function [13j . The parameter rg by its mean- 
ing is the strong interaction radius, usually taken as rg 1/m. Later on 
we cite it in all expressions, but as we checked, its impact on the results 
of fitting is small. Strictly speaking the expression ^ is valid in the region 
vmro < 1, however for considered kinematics it is a rather well approximation. 
To go under threshold one needs to perform the common replacement w — > z-D: 



~ 1 \/4m2 - M2 

log(2{}TOro) -h27 + 7rcot«) = (3) 

J Zm 

Due to asymptotic behavior ip{z) ~ log{z) this expression is finite at threshold 
Mc. Moreover it describes both all bound states (pionium atoms) and electro- 
magnetic interactions leading to unbound states under threshold [TU] . 



^More pessimistic view on impact of higher order terms has given in [7] 
^The photons radiation leads only to smearing of Coulomb contribution and will be treated 
elsewhere 

^The similar approach to — > Stt decays was developed by V.Gribov 2111] 
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Above the threshold the expression ^ acquires the imaginary part. Using 
the well known [13j relations for digamma function we obtain 



T = Re T + ilm t 
Rer = —a 



log(2mro) + 7 + ^ , c2\ 



= i_,-2.; (4) 

The above expressions allows one to calculate the electromagnetic corrections 
in wide kinematic region. For definiteness let us consider the fitting procedure 
in NA48/2. In order to extract the values of scattering lengths ao,a2 [2], the 
function 

F{M^)=N J f{ao,a2,m^)^M^,m^)dm^ (5) 

was fitted to the experimental mass spectrum. Here / is the theoretical dis- 
tribution from 5J, whereas <i> describes the experimental setup resolution and 
acceptance and obtained by Monte-Carlo simulation. Integral is replaced with 
a sum, and decay probability / is calculated in the middle of every M^^o^o bin 
with a width SM"^ = 0.00015(6*6^/0)^. The bins are centered in such a way, 
that the threshold ^rri^^ is placed exactly in the middle of one of them (here we 
will call this bin as the "central" one). 

This procedure is precise enough as long as one considers only strong interac- 
tions in final state. But the electromagnetic interactions lead to the sharp peak 
at the threshold, thus for the central bin the averaging of theoretical predictions 
must be done more carefully. 

So it is convenient to consider the behavior of the electromagnetic corrections 
in three separate regions. Up to the lower bound of central bin due to smallness 
of ^ in this region from ([3]) we get 

TTd TTCt [ G? 

T = —cot{ — ) + aUog{2dmro)+-f-^';{3) 
1 

This expression for r have to be used there instead of the velocity formula 

tor charged pions v = - — ^ . 

Above the charged pions production threshold, beginning from the top bound 
of central bin the interference between direct Tq and charge exchange T+ terms 
appears (unlike in [4]). As a result to take into account the electromagnetic 
effects in this region one has to add the interference term 

- 4aa, (^log{2mr„) + 7 + jr^^^ + ^^0.9(^2 + {^f)^ ToT+ (7) 

to the square of matrix element. 

For these regions the averaging procedure (just a value in the bin's center) 
exploited in [2l is precise enough. 
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As to the central bin, where the electromagnetic corrections come from two 
states of TT+TT" pairs (bound or unbound) the more accurate averaging is in 
use. Confining to the second order terms in scattering lengths, as it is done in 
previous consideration [S], substituting the relevant velocities p|4p . integrating 
the square of amplitude ([1]) in central bin and dividing result by the bin width 
<5M^, we obtain the mean value of correction to the square of matrix element in 
the central bin in the form 

= ^^!i(|lT2-4aa,(/o5(2vtmro) +7-0.5) T+To; 

Vt = v[Am- ^ ^ — ); uo = "0(4™ ) (8) 

The first term describes all bound states (pionium atoms). As to the second 
term it is a result of interference between the direct "unperturbed" production 
amplitude To and electromagnetic interactions leading to unbound states. Let 
us note that in ^ we cite only the main electromagnetic corrections in central 
bin. We estimated the possible contributions from unbound states (^ ) and 
interference between direct amplitude Tq and bound states {'^ T^Tq) in central 
bin. They are small (less then 1% in comparison with main terms in ([5])) and so 
can be safely neglected. 

Numerically the contribution of the first term (pionium atoms) from ^ 
into the central bin is near 2.5% of "unperturbed" decay probability, that 
coincides with the prediction evaluated from 14J. Simple estimates show that 
the interference of direct amplitude Tq with unbound Coulomb part (second 
term in ([S])) is significant and gives approximately the contribution of same 
order as bound states one. 

Now it becomes clear why the fit done in ^ with the pionium contribution 
as a free parameter leads to the Coulomb contribution, that is about two times 
larger, than can be provided by solely bound states. 

In conclusion let us stressed that at present we have the recipe allowing 
to introduce the Coulomb effects into the expressions for the square of decay 
amplitude (taken, for example, from [S]), averaged correctly over every bin, 
that can be multiplied by phase space factor, convoluted with the experimental 
acceptance matrix $ and compared with the experimental data. 
We are grateful to V. Kekelidze and J. Manjavidze who draw our attention to 
the problem and advice during all the work. 
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